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ABSTRACT. We prove optimal convergence estimates for eigenvalues and eigenvectors of 
a class of singular/stiff perturbed problems. Our profs are constructive in nature and use 
(elementary) techniques which are of current interest in computational Linear Algebra to 
obtain estimates even for eigenvalues which are in gaps of the essential spectrum. Further, 
we also identify a class of "regular" stiff perturbations with (provably) good asymptotic 
properties. The Arch Model from the theory of elasticity is presented as a prototype for 
this class of perturbations. We also show that we are able to study model problems which 
do not satisfy this regularity assumption by presenting a study of a Schroedinger operator 
with singular obstacle potential. 



1. Introduction 

In this paper we give sharp estimates for the asymptotic behavior of the spectral problem 
for the family of self-adjoint operators H K which are defined by positive definite quadratic 
forms 

(1.1) \) K (U,V) = t)b{u,v) + K 2 t) e {u,v), U,V G Q(f)b) C Q(f)e)- 

Here we have used Q(f)b) and Q(h e ) to denote the domain of definition of [){, and t) e and 
we assume that k 2 — > oo. Qualitative results for families of self-adjoint operators like H K 
have a long tradition. We are particularly influenced by the results from [28, 35]. Here by 
qualitative results we mean those results which prove (e.g.) that the spectral projections 
E K (-), H K = J A dE K {\) converge in some appropriate sense. 

To give a first idea of what is hidden within the abstract formulation ( 1 . 1 ) let us consider 
two simple examples that are representative for more complex model problems (studied 
later on in Section 5). The family of quadratic forms 

(1.2) f) K (u,v) = J u'v'dx + n 2 J u'v' dx, u, v € Hq [0, 2], n — > oo 
is paradigmatic for a regularly perturbed family, whereas the family 

(1.3) f) K (u,v) = J u'v'dx + n 2 J uv dx, u, v € Hq [0, 2], n — > oo 

is representative for the quadratic forms which violate our new regularity assumption. 
Note that in our relative theory the unbounded perturbation t) e in (1.2) is preferable to 
the bounded perturbation f) e in (1.3). Here we have used Hq(-) to denote the standard 
Sobolev spaces. 

Date: February 13, 2009. 

This work is based on a part of author's PhD thesis [14], which was written under the supervision of Prof. Dr. 
Kresimir Veselic, Hagen in partial fulfilment of the requirements for the degree Dr. rer. nat. The paper has been 
written at the time when the author was a member of "Institut fur reine und angewandte Mathematik", RWTH 
Aachen University. It has been submitted for a review since February 2008. 



2 



LUKA GRUBISIC 



The limit of the families like (1.2) and (1.3) can be a non-densely defined operator 
and we use the theory of [35] to study the convergence of such h K and associated H K as 
k — > oo. Let now the operator Hoc (in general non-densely defined) be the limit (in the 
sense of [35]) of H K as n — > oo . We use Af , i € N to denote the discrete eigenvalues 
of H K , which are below the infimum of the essential spectrum and are ordered in the 
ascending order according to multiplicity. By vf € Q(f) K ), H K uf = Xfvf and \\v" | = 1 
we denote accompanying eigenvectors. Here we allow n > or formally n — oo. Using 
the perturbation techniques from [13, 15, 16, 17] we prove (among other results) in the case 
of regular family of the type (1.1); for a definition see Section 1.2 below; the estimates 



and we compute the constants lb and ubj, i = 1,2,3 explicitly for several concrete model 
problems. Further, we also give a formula for determining a critical k such that (1.4)- 
(1.6) hold for k > n and we show that the estimates are optimal in the sense that 



To show that our abstract approach to problems (1.1) does not incur accuracy tradeoffs 
— when applied to concrete problems — we consider several case studies. A prototype for 
the (less trivial) regular problem is the Arch Model from e.g. [6, Chapter 8.8:3]. In our 
case study we compute explicit estimates for the asymptotic behavior of the eigenvalues 
and spectral projections of the low frequency problem as the diameter of the arch goes 
to zero. The limit of such family of arches is the so called Curved Rod Model from [19, 
32]. On the other hand, Schroedinger (like) operators from [2, 5, 9] are representative for 
(higher dimensional) operators which have less "well-behaved" spectral asymptotic. More 
to the point, in the case of the Schroedinger (like) operator from (1.3) we obtain the same 
optimality statements, but the convergence is of the fractional order 0(-^), a = | ( cf. 
[9, 13] for higher dimensional problems in unbounded domains). These concrete examples 
determine a framework for presenting our (otherwise) more abstract results. 

1.1. Local (resolvent) estimates. We approach this analysis by reformulating the conver- 
gence problem so that the perturbation framework and the error representation formulae 
(this is the main constructive feature of our framework) from [10, 13, 15, 16, 17] can be 
applied as a backbone of our construction. A difference between our approach and the 
standard results of works like [5, 7, 9, 26] can best be seen when considering a way to 
compute a constant ub 3 for an estimate like (1.6). The standard approach requires a study 
of the integral 



(1.6) 



(1.4) 



(1.5) 




D e K\ spec(Hoo), 





where C(A°°) is a circle in the resolvent set of H K which has A°° in its interior and the rest 
of the spectrum in its exterior. This frequently leads to cumbersome estimation formulae. 
Thanks to the local character of the error representation formula from [13], we are able to 
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base our theory on a study of the integrals 1 

1*00 

(1.8) K 00 ,H-\°°)-K°°,H- 1 0= / IIH^H-VH 2 rfr, i = l,...,m. 

Here H e is the operator defined by h e in the sense of Kato and m £ N is the multiplicity 
of A°°. The results from [3, 9] show that the integrals (1.8) are better amenable for a 
quantitative study than are (1.7). 

Due to the difficulties in dealing with a formula like (1.7), typical results from semi- 
classical analysis from e.g. [7, 28] establish only the fact that the projections converge in a 
much weaker sense (than is the convergence of spectral projections in norm) without giv- 
ing information on the speed of convergence as measured by the coupling k 2 . The nearest 
in spirit to our analysis is the approach of [26]. However, in this work only a particular 
family of model problems is considered and no estimates for the convergence of E K (-) in 
(unitary invariant) operator norm(s) are presented. Furthermore, the authors do not discuss 
the radius of convergence of their "asymptotic" expansions. For the geometric theory on 
the relationship between two projections and the importance of establishing convergence 
estimates for all unitary invariant operator norms we refer the reader to the seminal works 
[8, 18]. 

1 .2. A notion of regularity. Let us now make precise what we mean by the regularity of 
h e . In the terminology of [28] a family of the type (1.1) is said to be non-inhibited stiff if 
h e is a closed and positive quadratic form and the subspace 

(1.9) N(f) e ) := {u € Q(f)e) : M«] : = W"> u ) = °l 

(of H) is nontrivial. For technical convenience we assume ( without reducing the level of 
the generality) that f)& is positive definite and use and H e to denote the self-adjoint 
operators which are defined in the sense of Kato by and h e respectively. 

We identify the regular family of quadratic forms — with structure (1.1) — by requiring 
that hf, and h e satisfy a Ladyzhenskaya-Babuska-Brezzi type condition 

1 /2 

for some 6, t > 0. The condition (1.10) is equivalent with the claim that R(Hy 2 H b ~ 1/ ' 2 ), 

l<2 1/2 

the range of the operator H e ' H fe ' , is closed in H, cf. examples (1.2) and (1.3). The 
ramifications of the assumption (1.10) will enable us to formulate a new method for study- 
ing integrals (1.8) for this class of model problems and thus complement the study of 
singular obstacle potentials from [3, 9]. 

1.3. An outline of the paper. Let us finish the introduction by briefly outlining the struc- 
ture of the paper. In Section 2 we introduce the notation and present the qualitative conver- 
gence framework from [35] . The main approximation results of the paper appear in Section 
3. To be more precise in Section 3.1 we review the operator matrix approach to Ritz value 
estimation from [17, 13]. In Section 3.2 this approach to spectral estimation is specialized 
to the problems of the large coupling limit. In particular we make precise in which sense 
can these estimates be considered sharp. We also revisit, in Section 3.2.1, the example 
from [13] to show how do (1 .4)— (1 .6) look in praxis for a non-regular h e . In Section 4 



'The notation (■, •) and || • | always refer to the scalar product and the norm of the background Hilbert space 
H. The functions of the operator like H^ 2 are always meant in the sense of the spectral calculus. By PN(t) e ) 
we generically denote the H orthogonal projection onto the space N(h e ) 
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we characterize regular perturbations h e and give convergence estimates which utilize this 
additional structural information. In Section 5 we consider a model problem from the elas- 
ticity theory and show that its asymptotic behavior is regular. In the last section we put the 
results in the broader context and give an outlook of further research. 

At the end we would like to emphasize that our study is distinguished by its constructive 
character. This can be seen in the fact that we give a general method to compute the 
constants lb and ubj, i = 1, 2, 3 (as functions of H K and v°°) in (1.4)— (1.6). With such a 
result we give a method to establish both a first order correction for the limit eigenvalue 
A°°, as well as to assess the quality of this approximation to Xf. The optimality result is 
a justification of this claim. For other connections between the elementary linear algebra 
and spectral theory we refer the reader to [31]. 



2. Convergence of non-densely defined quadratic forms 

In this section we fix the notation and give background information on the previous 
results which we use. We follow the general notational conventions and the terminology 
of Kato [20, Chapters VI- VIII]. Minor differences are contained in the following list of 
notation and terminology. 

• Ti ... is an infinite dimensional Hilbert space, can be both real or complex 

• (•,•); || • || .. . the scalar product on Ti, linear in the second argument and anti-linear 
(when Ti is complex) in the first; the norm on Ti 

• Ti\<S)Ti2... the direct sum of the Hilbert spaces Ti\ and H2, for any x e Ti\ © Tii 

we have x = x\ © x 2 = for Xi <G Hi, i = 1, 2 

X2 

• spec(H), spec ess (H); A CSS (H) ... the spectrum and the essential spectrum of H; 
the infimum of the essential spectrum of H 

• A < B ... order relation between self-adjoint operators (matrices), is equivalent 
with the statement that B - A is positive 

• C(H); C(Tii 1 Ti2)... the space of bounded linear operators on Ti, which is equipped 
with the norm || • ||; the space of bounded linear operators from Ti\ to H2 

• RpO, N(X) ... the range and the null space of the linear operator X 

• ... the generalized inverse of the closed densely defined operator A. If A has 
the closed range then A^ = (A(A* A) -1 )* is bounded, see [24]. We will extend 
this notion below to hold for non-densely defined self-adjoint operators. 

• P, P±... the orthogonal projections P and P± := I — P 

• ... a permutation of N 

• diag(M, W) ... the block diagonal operator matrix with the operators M, W on 
its diagonal. The operators M, W can be both bounded and unbounded. The same 
notation is used to define the diagonal m x m matrix 

diag(ai, • • • , a m ), with Qi, • • • , a m on its diagonal. 

• s i(A) > S2(A) > ■ ■ ■ , s max (A), s m i n (^4) ... the singular values of the compact 
operator A ordered in the descending order according to multiplicity, the minimal 
(if it exists) and the maximal singular value of A 

• I X I ... a unitary invariant or operator cross norm of the operator X. Since 
I • I depends only on the singular values of the operator, we do not notation- 
aly distinguish between the instances of the norm ||| • ||| on C(Ti), £(R(P)), 
£(R(P), R(P)- 1 ), or such. For details see [30]. 

• tr(X), II X \\hs ■■■ the trace a the Hilbert-Schmidt norm of the operator X, it 
holds II X \\ HS = y/ti(X*X), see [30] 
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As a general policy to simplify the notation we shall always drop indices when there in no 
danger of confusion. 

Let us assume that we have a closed, symmetric and semibounded from below form h 
with the dense domain Q(h) C H as given in [20, (VI.1.5)-(VI.1.1 1), pp. 308-310]. The 
form () which has a strictly positive lower bound will be called positive-definite. This is 
also a small departure from the terminology of [20, Section VI. 2, pp. 310]. Such () defines 
the self-adjoint and positive definite operator H in the sense of [20, Theorem VI.2.23, pp. 
331]. Furthermore, the operator H is densely defined with the domain 2?(H) C 2(f)) and 
V{H^ 2 ) = Q(h). We also generically assume that H has discrete eigenvalues Ai(H) < 
• • • < A m (H) < • • • < Aoss (H) , where we count the eigenvalues according to multiplicity. 
Another departure from the terminology of Kato is that we use t)(ip, (f) to denote the value 
of () on ip, (f) e 2(f)), but we write h[i/>] := h(V>, tp) for the associated quadratic form h[-]. 
We also emphasize that we use •* to denote the adjoint both in the real as well as in the 
complex Hilbert space H as is customary in [20, Chapters VI- VIII]. 

In order to be able to handle the problems of the type (1.1), we shall need to work with 
operators that are not necessarily densely defined, cf. (1.2) and (1.3). We use the notion 
of the pseudo inverse of the operator H that is assumed to be self-adjoint in the closure of 
its domain of definition 2?(H) C ri (tacitly assumed to be a non-trivial subspace). A 
definition from [35] will be used. The pseudo inverse of the operator H is the self-adjoint 
operator H defined by 

©(H) = R(H) © D(H)- 1 -, 
U(u + v) =H-y u e R(H), v e 25(H)- 1 -. 

It follows that H = H 1 in R(H) "'" and H is bounded if and only if R(H) is closed in H. 
When considered solely in 2?(H) " " the operator H is obviously self-adjoint, so we can 
also use the spectral calculus from [29] to define the generalized inverse, which extends 
the definition from the case of the densely defined operator, as 

V(H^) = {ueH : J f 2 (\)d(E(\)u,u) < oo}, 

where E(-) = E H (-)P V{H) . Obviously, we have 2?(H) N (H) = D(H*) and the identity 
H^m = Hw, u e ^(H 1 / 2 ) holds. In further text we shall tacitly drop the notational 
distinction between the generalized and pseudo inverse. The usual monotonicity properties 
can be extended to the generalized inverse. In particular it holds 

(2.1) IIh} 7 2 u\\ < \\Ul /2 u\\, u £ P(H2 /2 ) & ||H2 /2t w|| < ||H} /2t w||, u e V(U\ /2 ). 

This monotonicity principle is the main ingredient of the proof of the convergence result 
for (1.1). When dealing with non-densely defined forms this principle can be formulated 
as follows. Let hi and f)2 be two closed positive definite forms and let Hi and H2 be the 
self-adjoint operators defined by hi and h 2 in Q(f)i) and <2(f) 2 ). We say hi < h 2 when 
Q(h 2 ) C Q(hi) and 

(2.2) hiM = ||h; /2 u || 2 < h 2 M = HH^H 2 , u g Q(h 2 ). 

Equivalently, we write Hi < H 2 when hi < h 2 . Now, we can write the fact (2.1) as 

(2.3) Hi < H 2 <^=> H 2 < H\. 
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Let us define, for non-inhibited (see definition (1.9)) quadratic forms like h K from (1.1), 
the domain Qoo := {u e Q : Hindoo h K [u] < oo}, then according to [29, 35] the 
symmetric form 

f>oo(w, v) = lim t) K (u,v), ti,wefioo 

K — »00 

is closed in and it defines the self-adjoint operator Hqo there. Further, it holds that 

= s-linifc^oo H~ 1 . The general framework for a description of families of converging 
positive definite forms will be the following theorem from [35]. 

Theorem 2.1. Let s n , h„, u„ and hoc be closed symmetric forms in H such that they are 
all uniformly 2 positive definite. 

(1) Ifs n > s„+i > hoc where 

t)oo(u,v) = lim s n (u,v), u,v E I J Q(s n ) 

then hoo is closed with Q(hoo) = [J„eN 2( s «) "°° and H L = s-lim n 

(2) Ifu n < u n+ i < where 

t)oo(u,v)= lim u n (u,v), «,»eQ(l)oo) 

n— >oo 

then is closed with Q(f)oo) = {/ € f] n efi 2( u «) : supu n [/] < oo} and 
Hj = s-limn Ut. 

(3) If u„ and S n are as before andu n < h„ < s„ a/so holds, then 

i)oo(u,v)= lim h„(u,w), Q(f)oo), 

n— >oo 

Hj,, = s-limHt . 

For the families of forms which satisfy the assumptions of Theorem 2. 1 the following 
qualitative convergence result on spectral families has been established in [35]. 

Theorem 2.2. Let h„ be a sequence of positive definite forms that satisfies any of the 
assumptions of Theorem 2.1 for u n , S n or h„. Let there also be the positive definite form s 
such that h„ > s and A e (S) > 0. Then 

(2.4) \\E n (D) - E^WW ^0, D < A e (S),D ^ spec^). 

The results like Theorem 2.1 have independently been obtained in [28, 29]. We have 
opted for Theorem 2.1 since it extensively uses the monotonicity (or "sandwiched" mono- 
tonicity) to establish the stability of the converging eigenvalues and this fits neatly into 
the perturbation framework of [10]. This was the chief source of motivation for the main 
construction from the PhD thesis [14] (those results appeared later in [13, 15, 16, 17]). 

3. A CONSTRUCTIVE APPROACH TO ASYMPTOTIC EIGENVALUE/EIGENVECTOR 

ESTIMATES 

Let us reiterate that we use the notion of the constructiveness in this paper in two con- 
texts. First, it should emphasize that all of our theory is bases on the error representation 
result like (3.8)-(3.9), below. But second, it is also meant to emphasize that in a result 
like those of the type (1.4)— (1.6) we present a way to construct an improvement to the ap- 
proximation A°° (of the eigenvalue Af ). The constants lb and ub,, i = 1, 2, 3 are explicit 
functions of the approximation defects ^(P), to be defined below and it is the aim of this 
section to reveal this dependence. 



'By this we mean that they have a uniform positive lower bound. 
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3.1. Background information on the block-diagonal part of the operator/form. In this 
section we review the results form our previous work which we use to prove our first 
contribution in Section 3.2. A reader who would like to go straight to the new results can 
do that directly after reading equation (3.1) and Definition 3.1 below. 

In this section we assume that we have a fixed closed symmetric and densely defined 
form h. We will review the basic spectral properties of the block-diagonal part of h with re- 
spect to orthogonal projection P, R(P) C Q(h) as is presented in [15]. In order to simplify 
the presentation we temporarily suppress (in the notation) the dependence of quantities on 
H, where there is no danger of confusion. Assuming that R(P) is finite dimensional we 
define the block-diagonal part of f) by setting 

(3.1) t) P (u,v) :=t)(Pu,Pv) + i,(P ± u,P ± v), u,ve Q(h P ) := Q(h). 

Obviously the form hp is closed and positive definite and so it defines the self-adjoint 
operatorHp in the sense ofKato. We further have (for a proof see [13, 15]): 

(3.2) R(H _1 -Hp 1 ) is finite dimensional. 

(3.3) W (P) := sup IfrM ~ < L 

ueQ(fi) "JpM 

A first consequence of these two features is the stability of essential spectra, namely Weyl's 
theorem gives spec css (H) = spcc css (Hp). Further, we have the estimate — of the same 
form as (3.3) — for the eigenvalues A,(H P ) and Aj(H), i e N which are below the infimum 
of the essential spectrum A CSS (H) = A css (Hp) 

(3.4, ^-W <un i£N . 

The attractiveness of interpreting the form h as a perturbation of its block-diagonal part 
lies in the fact that 

(3.5) spcc(Hp) = spcc(S) U spcc(W) 

where S = (H 1 / 2 P)*(H 1 / 2 P)| R( . p ^ is a finite dimensional operator and W is the self- 
adjoint operator which is defined in R(Pl) by the quadratic form f)(P_L-, Pj_-). Since 
spec(S) is computable, we can start building our constructive estimation procedure on this 
fact. As a convention we will use fii < • • • < MdimR(p) to denote the eigenvalues of S. 
The numbers \ii will be called the Ritz values from the subspace R(P). In this section we 
also use the notation Aj := Aj(H). 

Let us now assume that dim R(P) = m £ N. To examine the relationship between h 
and hp in further detail define 

(3.6) Vi (P):=[ ^ax^ mm { | ^ € S, ||V|| = l}j , 

dim (5)— m— i+1 

for i = 1, . . . , m. It has also been shown in [15] that r/ max (P) = r] m (P). Although the 
perturbation <5p(h) := h — hp is in general — forsomeP, R(P) C Q(h) — not representable 
by an operator, the quadratic form <5 P (h)[-] := l][H p 1 '' 2 -] — hpfHp 1 ^ 2 -] can always be 
represented by the bounded operator block-matrix (with respect to P © P^ =7) 

V)' and (;6p(H)-) = 6 P (t,)[-]- 

Furthermore, [13, Lemma 2.1] gives that Si(T) = rn(P), i = 1, . . . ,m. The analysis 
of [13] now yields the conclusion that the test space R(P) can be used to generate good 
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approximation for the eigenvalues Aj, i = q, . . . , q + m — 1 when f] m (P) is smaller than 
half of the relative gap 

( \+m — A*m Ml — \-l \ 

Iq ■= mm i v i ' ~\ — f' 

\+m + Mm Ml + J 

Ample numerical evidence corroborate that such estimates are robust (with regard to scal- 
ing) and sharp. Assume that t] max (P) < \^ q and that dimR(P) = m, where m is the 
multiplicity of the eigenvalue X q . Using [13, Theorem 3.3] we conclude that the operator 
matrix 



(3.7) S' P {H q ) 



I - AqS" 1 r* 

r i-XqW- 1 



which is the block-matrix representation (with respect to P © P± = I) of the quadratic 
form 

S s P (i) q )[] ■= h(H p 1/2 -,H p 1/2 .) - A 9 (H p 1/2 -,H p 1/2 .), 

satisfies dim N(S s p (H q j) = m and the mechanism of [31, (1.1)— (1.2)] — also known in the 
Linear Algebra as the Wilkinson's Schur complement trick (see [27, pp. 183] and [13, 
Theorem 3.3]) — allows us to conclude 

(3.8) I - A«jS _1 = r*(i — A«jW _1 ) _1 r 

(3.9) = pt + x q r*w-^ 2 (i - A 9 w- 1 )- 1 w- 1 / 2 r. 

Identity (3.8) is the basis of the proof of [13, Theorem 3.3] which we now quote. Note 
that (3.8)-(3.9) also hold for X q which is in a gap of the essential spectrum. Based on the 
definition (3.6) we now define (for later usage) the approximation-defects for h K . 

Definition 3.1. Let the sequence h K be given and let the orthogonal projection P be such 
that R(P) C Q(h K ) and dim R(P) < oo. We write P) for 7]i(P) from (3.6) when 
applied on h K . We call ^(k, P) the k- approximation defects. If we are given a subspace 
^3 = R(-P), then we abuse (simplify) the notation and freely write r]i(K, *}3) = /^(k, P). 

Theorem 3.2. Let the discrete eigenvalues of the positive definite operator H be so ordered 
that X q -\ < X q = A«j +m _i < X q+m . Let R(P) C Q{h) be the test subspace such that 
dim R(P) = m and 1 ^^ P ) < 7<j- Then we have 

(3.10) I diagt^^)™! I < I diagMP))^ I . 

Mi 9q,r] m (P) 

where Q q , c := max { c ' c , c _ c ) for q > 1 and we set 

(1+ 1_£ Mg— 1 \1 )*q + m 

01, C := 01 := A m +i+^ m " He re we use diag(ai)™ 1 to denote the m x m diagonal matrix 
with scalars ai on its diagonal and ||| • ||| denotes any unitary invariant matrix norm and 
Hi are the Ritz values from R(P). 

In the case in which we do not have explicit information on the multiplicity of X q 
we have a weaker upper estimate. There is also an accompanying lower estimate which 
establishes the equivalence of the estimators r)i(P) and the error. Assuming that H = 
/ A dE(X) and that we use Vi, Hvi = XiVi, \\vi\\ = 1 to denote eigenvectors and ipi £ 
R(P), Etpi = [iiipi, = 1 to denote Ritz vectors, we collect some representative 

spectral estimates (bases on R(P)) from [13, 17]. 
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Theorem 3.3. Let the discrete eigenvalues of the positive definite operator H be so ordered 
that X m < X m +i and let A Sl < \ S2 < ■ ■ ■ < X Sp be all the elements 3 o/spec(H) \ {A G 
spec(H) : A > A m+ i}. If jj^ypj < X \ +1 j rl ^ m then there exist eigenvectors Vi, 
Hv-i = XiVi, \\vi\\ = 1 and Ritz vectors ipi € R(-P), B.ipi = ^iipi, \\ipi\\ = 1 such that 

(3.11) 1 E(p, m ) - P I < , 

A m +1 - Mm V 1 ^ Vrn{P) 

m m I \ I 1 m 

(3.13) |k-^|l< max ^ 



Ae S pcc(H)\{A s } \\-m\ y/l-ThniPy 

(3.14) — -j— — = \\Vi-ipi\\ H , i = l,...,m. 

Wi\ Xi 

Here P Si is the orthogonal projection onto the linear span of {tjjj : j = YL\=i m k + 
!)•••! y^l— i nik} and m, is the multiplicity of the eigenvalue X Si , i = 1, . . . ,p. Obviously 
the identity P Sl © P S2 © • • • © P Sp = P holds. In the case in which Ai = A m we can drop 
the constant j^— from the lower estimate. We can also allow for other cross norms ||| • ||| 
of the diagonal matrix diag((7?j(P)) in (3.12). 

The proof of the estimate for the spectral projection (3.11) can be found in [17], the 
proof of (3.13) is in [13] and identity (3.14) is well4mown. For reader's convenience let 
us also point out that the problem of estimating the spectral projections E(J) — where 
3 is some contiguous interval whose boundary points are not the accumulation points 
of spec(H) — can be seen as problem in obtaining a robust computable estimate of the 
Cauchy integral 

(3.15) \\E(3)-P\\ = ±\\<f (C-HpV-MC-H)- 1 ^. 

By €(3) we denote the circle in the resolvent set of H such that 3 is in the interior of the 
associated disc and the rest of the spectrum is outside the disc. However, contrary to the 
intuition, the direct analysis of (3.15) is not the most natural way to obtain computable and 
robust estimates of ||-E(3) — P||. A problem is that, although the integral of the resolvent 
difference does not depend on the integration path €(3), estimates of it do. Furthermore, 
the circle is only one of many possible curves which should be taken into account. As 
an alternative we consider the approach of the (weakly formulated) operator equations. 
Not only are the estimation formula which are so obtained sharp (see [17, Remark 2.3]), 
but also the technique allows for a natural consideration of estimates which utilize other 
operator cross norms ||| • |||. Such results are known as sin 9 theorems in the recognition 
of the milestone work [8] and have been extensively studied in the computational Linear 
Algebra, see [22, 21] and the references there. We use a recent generalization of those 
results, which is particularly suitable for an application in the quadratic form setting, see 
[17]. 



'We assume that 1 < si < S2 < • • • < s P < m. 
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Remark 3.4. Note that as r] Si (Pi) — > we have Q Si ,rj m (p,) - * mm ' T~ i T7 ^ ' A 'a ^ L l 
and min{g s . r)m . (p.) : i = 1, . . . ,p} quantifies the minimal relative gap among the eigen- 
values A S1 < A S2 < ■ ■ ■ < X Sp . Note that the relative gap Sij ^ s .(p 4 ) distinguishes better 
between the close eigenvalues than the absolute gap, eg. min{A Si+1 — X Si , X Si — X Si -i} 
is an example of an absolute gap. In Theorem 3.3, equivalently as in [11, Proposition 2.3], 

we have that when r) mi (Pi) < k mirife^- ^ 3fc i = l,...,p then 

1 < 3 



3.2. Estimates for the spectral asymptotic. We will now use Theorem 3.3 to obtain con- 
vergence rate estimates for (2.4). This is the central result which guaranties the stability 
of the spectrum of the converging family of forms fj K . Subsequently we will also prove 
results like (1.4)— (1.6) and use the motivating example of the Schroedinger operator with 
a singular obstacle potential from [13, Section 4] to show our estimates in action. 

Although we are working under the assumptions of Theorem 2. 1 , we assume — in order 
to be more explicit — that we have the non-inhibited stiff family h K from (1.1). The form 
hoo obviously defines the self-adjoint operator Hoo in IM(h e ). By Hoc = J XE CC (X) we 
denote the spectral representation of Hoo in N(h e ). We identify S (X) (-) with £ 00 (-)P|\|([ )e ) 
and write Hoo = / XE OQ (X) for the non-densely defined — in the space H — operator Hoc. 
Let 3 be a contiguous interval in M., then (S^ := R(-EooP)) is a subspace of Q := Q(f)b). 
Let now 3 be such that (E^, is finite dimensional, then K-approximation defect is given by 

(3.16) »>*(«,<!*):= maxmin f' °° V> g 5, M = l} , 

dim (oj— m— 2+1 

where H~3 := (Hj^^wj) = (H~ 1 ) Eoo ( 3 \ and z = 1, . . . , dim R(€^). To further 
simplify the notation we set 3) := ^(k, <B^ x> ). Theorem 2.1 now obviously yields 
lim r\i (k,3) = 0, i = 1, . . . , dim R(€^). 

K — »00 

Similar construction can be performed in the case in which (E^, is infinite dimensional. 
The main features which are lost in this generalization are the easy computability of 
spec(S^T 1 ) — spcc(Hg3 Eoo(3)), the property that always r] max (K, 3) < 1 and the re- 
sult on the stability of the essential spectrum. This makes, in general, such method less 
attractive for practical constructive considerations. 

Let us first give a quantitative version of Theorem 2. 1 which is based on the application 
of Theorem 3.3. As a notational convenience we use Xf < ••• < X°° < X^ s and 
XI < ■ ■ ■ < Xf < A£ ss to denote the discrete eigenvalues below the infimum of the 
essential spectrum of the operators Hoo and H K respectively. 

Theorem 3.5. Let H K = J X d E K (X) be the operators which are associated with the 
family of forms fj K . Take D G R such that X^ < D < A^ +1 and set 3 = (— oo, D], then 

(3.17) J7i(/c,J)<l, i = l,...,m, 

(3.18) j j <v m (K,3), j = l,...,m, 



(3.19) \\E K {D) - E^D)]] < ' / "' (,V ' J) 



\D-X%\ v/1-7^(k,J) 
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for k large enough. (For the meaning of the phrase large enough see Remark 3.6.) 

Proof. Statement (3.17) is a direct consequence of [13, Lemma 2.1]. Let us now remember 
(3.4). This estimate is the consequence of [15, Theorem 4.5] which, when applied to the 
form h K and its €oo block-diagonal part (f^)^, yields 

(1 -?7 m (K,a))(h K ) Coo < h K < (1 +?7 m (K, J))(h K ) Coo . 

Let (H K ) c be the self-adjoint operators which represent the forms (f) k ) g in the sense of 
Kato,thenAf e spec (H„) ^. SetW K = (H K ) <Sgo E 00 (D)± and = H 00 E 00 (D)± 
then Theorem 2.1 implies that = s-lim K ^oo Wj,. By the construction of W K we 

have R(S 00 (D)) _L w for any w e V(W K ). This implies Ai(W«) -» Ai(Woo) = A~ +1 . 
On the other hand, since 

spec((H re ) ) = {A e spcc(Hoo) : A<L>}UspecW K 

it follows that there is k suc h that 

[\™,D] cR\spcc((H K ) £ J, k>« . 

Since T) m (n,3) — > we conclude that for k > k (here we slightly abuse the notation) 
the estimate iq m (n,3) < ^ D+A S holds. Now, the conclusion (3.18) follows from [15, 
Theorem 5.2]. Equivalently, the conclusion (3.19) follows from (3.11) and [17, Theorem 
3.2]. □ 

Remark 3.6. The coupling constant ko is large enough when 

1 \°° — \°° 
Vm ( K ,3) < 1 Am+1 Am 



J A rn+1 "+~ A m 

for k > k . This follows by a similar consideration as in Remark 3.4. 

A direct application of the results from [17, Section 3] and the results of Theorem 3.5 
is the following corollary. 

Corollary 3.7. Assuming the setting and the notation of the previous theorem we have 

^DX% I diag(( % (K, J))^) 8 diag((Tfc(«, J))^) 1 



-£«,(£>) |||< 



|£>-A° 



V 1 - Vm(K,3) 



In the case in which 3 = [D-,D+] and X^_ 1 < D- < A£ < A^ +m _ 1 < D+ < A£ +TO , 
k> kq then 



(3.20) 11^(0) - £7ooP)|| < 



|£>+-A~| |Af-D_|J ^/l-^J) 



An easy comparison with the single operator estimates from Theorem 3.3 reveals that, 
unlike the spectral family estimate (3.19), the eigenvalue result (3.18) is suboptimal in the 
asymptotic setting. The problem is that we can not uniformly apply the estimate (3.12) 
on all the operators H K , k > k since we have no information of the multiplicity of the 
eigenvalue Af for all k > k - We only know the multiplicity of A?°. The only statement 
which we can make in general is a lower estimate on the convergence rate. A way to 
solve this multiplicity problem will be presented in Section 3.2.4, for now we only give the 
following result. 
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Corollary 3.8. Assuming the setting and the notation of Theorem 3.5 we have 

\ OO m m |\K \OOl 

2A ™ i= i i= i A * 

Furthermore, for each k > Owe can chose eigenvectors vf, H K wf = Xfvf, \\vf\\ = 1 and 
v°°, Hoo-yf 5 = Xfvf, \\vf\\ = 1 such that 



E^ 2 m<E 



2A~ " tt ^^°°] 

One situation in which we can readily obtain upper estimates like those from Theorem 
3.3 is the case when we know that A°° has the multiplicity one. This is frequently a case for 
the ID differential operators. Also, the lowest eigenvalue of many Schroedinger operators, 
like those from [7] have multiplicity one. In what follows we use || • || A -i = || A -1 / 2 • | 
to denote the standard A _1 -norm, which is associated to a positive definite operator A. 

Theorem 3.9. Assume the setting and the notation of Theorem 3.5, and let A^°, q e N be 
of multiplicity one then 

(3.21) lim 2 /\ . - 1, 

W<-<°] 



(3.22) lim 2 / ;L = 1 

k^oo r]i(K, A^°) 

Proof. By the same argument as above we may assume that we have kq such that 

-i \oo \oo \oo \oo 

Vi (k,X™)< ^minl X 1 7 «> Ko . 

/ J- v : q f — q L \oo i \ oo 7 \ oo i \ oo J ' u 

Theorem 3.5 yields that there exist D-,D + such that < D- < A^° < D + and 

(3.23) X*_ 1 < D- < Xq < D + < A£ +1 , k > k - 

According to [13] we conclude that we may apply the error representation formula (3.9) to 
the operator H K and the test vector such that HocV™ = X^v^, \\v^°\\ = 1. To the 
vector we can define the residuum as the functional t£ := H K u^° — A^°u^° and the 
identity 

H^Hk)^ =K CO ,H co ^) ?? 2 ( K ,A-) 

can be established by an easy computation. Also note the following identities 

„ M| |h K K^)-(h K ) goo K^)| 
C m N-i = max H : ' — = max ----- - — ------ 



' M t H -)eoo »£S\{0} IKH^J/fvH "£S\{o} ll(H K )g5 

IM«>«S )-(Weo>>''2 )l 



max 



"£fi\{o} ||(H K )l /2 w|| 
i>_LN(f, e ) llv k;c =° 11 

where (•, •} is the standard duality product. Analogous manipulation and the error repre- 
sentation formula (3.9) yield the conclusion 

Si -i, ^ ((H K ) e j^(I-A-(H K )^)-HH K )^r«) 
t? 2 (k,A~) A- »^(«,A~) 

= 1 + 0(||(H K )^ /2 P N(W J| 2 ). 
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Finally, Theorem 2.1 implies (3.21). The conclusion (3.22) follows from (3.14) □ 

We would like to emphasize that in this result the monotonicity of the family h K played 
a role. It is possible to prove the result without the property A^° > A£. The proof is 
technically more involved and it does not further the understanding of the problem, so 
we leave it out. This theorem establishes that the estimate — which follows directly from 
Theorem 3.2 — is sharp. We formulate this as the following corollary. 

Corollary 3.10. Assume the setting of the preceding theorem then 



A„ m n-f q+ 2_ _2 til 

This estimate is sharp in the sense of (3.21). 

3.2.1. A concrete example. Let H K be the operators which are defined by the family of 
positive definite forms 

/•oo />oo 

(3.24) f) K (u, v) = / d x ud x v dx + k 2 / uv dx, u, v e Hq (K+). 
Theorem 2.1 readily yields 

f)oo(u, u) = / d x ud x v dx, u,v € Hq[0,1]. 
Jo 

Here we have used Hq[0, 1] and (M+), M+ := [0, oo) to denote the standard Sobolev 
spaces. We also identify the functions from Hq[0, 1] with their extension by zero to the 
whole of K+ and write H&[0, 1] C #o(M+). We a l so formally write H K = -d xx + 

« 2 X[i,oo) and Hoc = -d xx and chose 

, x [V2sin(fc7ra;), 0<x<l 

(3.25) Ui (x) = r v ~ ~ 

10, 1 < x 

as a test function(s). A simple computation yields that Xf is a solution of the equation 



(3.26) Vk 2 - X K = - cot(V>) 

and we know that each Xf has the multiplicity one. The quotient Al A ~ Al can be represented 
(for k — > oo) by a convergent Taylor series (see [25]) 

A ?°- A ? -O 1 o 1 , of 1 , 1 2A 1 ln /l J 4 A 1 



(3 ' 27) ^- 2 «- 3 ^ + 8 U + 4!^^- 10 U + i!^^ + ---- 

Using the Green functions we also directly compute i] 2 (k, X°° ) := . For computational 
details see [14]. 

By utilizing the information from (3.26) we can establish 



(3.28) (i_^_l_)47r 2 =:^( K )<A 2 (H), k > 5, 
which leads, in combination with (3.12), to the estimate 

(3 .29) _^<^<£M±i!_^ = i£ + 'o(i), «> 5 . 

Similar sharp results can be obtained for other A°° and using (3.14) for corresponding 
eigenvectors. We tacitly leave out the details. 
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3.2.2. A remark on higher dimensional singular obstacle problems. This paradigm has 
been applied in [14] to operators which are defined both in H 1 (M n ) as well as in 
where Cl C R n is a bounded domain. The only ingredient which is necessary is a result on 
the behavior of the momenta 



(3.30) 



(* OO 

(/, H" 1 /) - (/, nlf) = / IIH^H-VH 2 dr, f £ Co 

Jk 2 



Estimates of such momenta have been obtained on many places in the literature. We illus- 
trate our point by a consideration of a model problem of the electro-magnetic waveguide 
O x WL, where the section O C R 2 is a smooth and connected domain. The material 
C O of very large conductivity is compactly immersed in O, which is to say that the 
closure cl(f2) is contained in O and that Cl is bounded. The dielectric material is now mod- 
eled by U = O \ Cl. Assuming that the boundary of Cl is sufficiently smooth we study the 
eigenvalue problem for H K = —A + K 2 xn- Here, \n is tne characteristic function of CI 
and H K is the operator which is defined in the sense of Kato by the quadratic form 



f • + n 2 f X Mp<t>, V, e Qoo := Hi (O) 
Jo Jo 



where k £ M + and V is the usual gradient operator on Hq(O), the Sobolev space of 
functions with zero trace on the boundary dO. 

This problem has been analyzed in [12]. Let us assume that < D < A~ +1 , for 
some to £ N. We compare P = E^X™) and Q K = E K (X%), where H K = j X dE K (X) 
is the spectral integral in L 2 (0) and = J X dE OQ (X) is the spectral integral in L 2 (U). 
Since, as has been shown in [12], 

(v°° H^ 1 v°°) - (v°° H f v°°) = - — - — / 

l A i ) Jd 

we have coarse eigenvector estimates 

^- p ^j^%T K + °^^^^ 



i i r dv? dv? , i , 
dn dv dv + y ^i 2 ' 



(3.32) min ^LBI J_ < < 4 1 



K 



which can be improved in a straight forward manner by bringing the factor J dn 
into estimates, as has been shown in [12, Section 2.1]. The last inequality in (3.31) and 
(3.32) hold for k large enough. The optimal eigenvalue estimate can easily be constructed 
from Theorem 3.3 and 3.9 and we know that the eigenvector estimate (3.32) is optimal in 
the sense of (3.14) and (3.22). Remark 3.6 indicates how to assess the radius of conver- 
gence of these first order estimate(s). 

3.2.3. Remarks on (finite) eigenvalues in gaps of essential spectrum and on general con- 
verging families h K . We have said that the theory can be applied to eigenvalues which 
are in the gaps of essential spectrum. Since we do not consider any model examples which 
show such behavior (e.g. operators with periodic boundary conditions) we will only briefly 
outline a possibility to obtain results like Theorem 3.9 or Theorem 3.5 in this setting. 

In dealing with the eigenvalues in gaps of the essential spectrum we do not have the 
safe convergence environment of Theorem 2.1. Instead, we have to have an a priori infor- 
mation that the assumption like (3.23) holds. An example of how to obtain this type of a 
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priori information can be seen on the proof of (3.28). To this end we would like to empha- 
size that such type of "precise" result on the separation of the target eigenvalue from the 
unwanted component of the spectrum is an unavoidable ingredient of all constructive spec- 
tral estimates. An assumption like (3.23) is equivalent to requiring that eigenvalue X q be 
stable under the perturbation h K , see [20, chapter VIII.4, pp. 437]. For a characterization 
of perturbations for which this holds see [23] and references therein. 

Given such an estimate — i.e. assuming that \ q is a stable eigenvalue — the appropriate 
result from [13] or [17] can be applied to obtain convergence estimates. We also empha- 
size that the theory of [13, 17] allows for more general spectral intervals 3. To be more 
precise, to establish an estimate like (3.20) the spectrum in 3 n spec (Hoc) does not have 
to be discrete. However, in such situation we have no guarantee that ?7max(K, 3) < 1 and 
obtaining computational formulae requires much more technical work. The precise use in 
a given situation is application dependent, but always follows the procedure outlined in 
Theorems 3.5 and 3.9. 

In the case in which we consider a general converging family of quadratic forms from 
[35] we cannot conclude that (H^J^p) = (H^ 1 so we have to use explicitly 

computable S^ 1 := (H K )~^^\ E .„ in (3.16) instead. If we set /xf := Aj(S re ), then 
substitutes for A°° in eigenvalue estimates like (3.18), (3.21), whereas the estimates for the 
spectral projections like (3.20) remain unchanged, e.g. we have the convergence estimate 

\\0 nil < \/ A m+l/C Vm (K,P) 



3.2 .4. A method to solve the multiplicity problem. A tacit assumption in this semiclassical 
analysis is that the operator Hoc is a well known object. In order to be able to apply 
Theorem 3.3 one should establish that there exists k > such that 

(3.33) Xg_ 1 < D- < Xg = X q+m _ 1 < D + < Xg +m , 

for k > Ko- However, if to > 1 it is not plausible to expect that (3.33) will hold in general. 
Instead, we will get a tight cluster of to eigenvalues (counting the eigenvalues according to 
their multiplicity) that converge to A^°. Since we aim to express the spectral information 
about H K in terms of the spectrum of Hoc we further opt to give specific values for _D_ 
and D + as functions of the gaps in the spectrum of Hoc. 

Theorem 3.11. Let the eigenvalues of the operator Hoc be so ordered that X q < ^ 1 < A^° = 
A g+ro-i < ^q+m- Define the measure of the relative separation o/A^° from the rest of the 
spectrum o/Hoo as the number 



\oo _ \oo \oo \oc 



I \oc i\oo'\oo_i_\oo I 
[ A q+rn + A q A q + A q-l J 

There exists k > such that for n > k 

3^0^°°) 

\ K \°°\ —, q 

V-M) <l] m (K,A q ) 3 ( A oc), ,m. 

1 7.(A~) 
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Proof. Since i] m (k, A^°) — ► 0, an argument analogous to the argument that led to Theorem 
3.5 implies that we can pick kq > such that for n > «o 

(3.35) 7? m («,A~) < i 7s (A^°) 



(3.36) |A^-A^| < - 7s (A^)A^°, fc = g, <? + 1, . . . , q + m - 1, 



1 

3 

(3.37) |C-A fc (H„)| > i 7s (A~)A fe (H K ), fc {«,..., g + m - 1}, C € C(A~). 

Here £(A^°) is the circle in the complex plane with the radius §7s(A^°)A^° and the center 
A^°. Assume k > k is fixed, then define the family 

(3.38) o(r) = (t) K )p + rS P (t, K ), r e C. 

This is a holomorphic family of type (B) (for the definition see [20, Chapter VII]). We know 
that 

(3.39) \6 P (i, K )[u]\<ri m {K,X^)(l) K )p[u}, u e Q, 
so [20, Theorem VII-4.9 and (VII-4.45)] imply that the resolvent 

i?(r,C) = (A(r)-CI)- 1 

can be represented by a convergent power series in t for ( E £(A£°). The power series for 
i?(r, C) converges for every 

(3.40) \r\<r = -^-^ inf JA_^ = _L^I 7s ( A °°). 

AGspoc((H K ) P ) 

In particular, assumption (3.35) implies that the series converges for r = 1. 
Define 



%):=-J-A(r)/ fl(r,C)dC, 



|A i (r)-A~|<- 7s (A~)A~, |r|<r . 



C(A~) 

then -B(t) is a holomorphic operator family and there exist m holomorfic functions Ai(r) 
such that Ai(t), • • • , A m (r) are all the nonzero eigenvalues of the operator B{t). Due to 
the assumptions we have made it follows that for i = 1, . . . , m 

1 

3 

Cauchy's integral inequality 4 for the coefficients of the Taylor expansion implies, for every 
i = 1, . . . , m, the estimate 

where A;(t) = A?° + r\ (1) + r 2 Af } + r 3 Af } + • • • . This yields 



\%(r) - A~ - rA« I < H < 



T 



1 



ro r -\t\ rl \ _ ill 

ro 



4 For further details see [1, Section 8.1.4] and [20, Section H~-3]. 
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for \t\ < ro- In particular for r = 1 there exists a permutation such that Aj 4 (l) = 
AJ +i _i, i= l,...,mso 

|\k \oo TWl ^ /• \oo^\oo ^Vm{ K t ) 1 

|A 9+^-l ~ \ ~ A ii I < \ >\ 7c (\oo) ~ _ 3r, m (« A~) • 

With this is the proof of the theorem finished. To see this note that it was established, in 
[20, (VII-4.50)], that A^ are the eigenvalues of the matrix Mk p = <5p(h K )(ufc, u p ), where 
Ufe, k = 1, . . . , to form an orthonormal basis for R(Eoo[D-D+]). Since 

( 5p(h K )( u ,«) = h K (P ±M ,P7i) + h K (P u ,P ±M ) = 0, R(P), 

we obtain A^ = 0, i = 1, . . . , m and the conclusion follows. □ 

Remark 3.12. The estimate of this theorem is optimal in the sense of Corollary 3.8. The 
upper estimate which has a similar form to (3.12) can be established for the limit eigen- 
values AJ° < • • • < A^. The the role of the constant from (3.12) is taken by the constant 
7min(A^) := min{7 c (A°°) : i = 1,.. ., to}, as is given by the repeated application of 
Theorem 3.11. We leave out the technical details. 

4. Spectral asymptotic in the regular case 
We now concentrate on the non-inhibited families 

(4.1) i) K (u,v) = t) b {u,v) + n 2 t) e (u,v), u,vGQ:= Q(h 6 ) C Q{t) e ), 

1/2 —1/2 

which satisfy the additional regularity assumption that the range of the operator H e H 6 
is closed in Ti. As already mentioned in Section 1 .2 this is equivalent with 

(4.2) || (Hi/ 2 H- 1/2 )t || =t<oo. 

With this additional requirement, which has a flavor of Linear Algebra, we can use an 
adaptation of the Lagrange-Multiplier technique to establish an upper estimate for the mo- 
menta 

/•oo 

(4.3) (f,H-\f) - (f,Ulf) = / IIH^H-VH 2 dr, f e Qoo := fi(M- 

The lower estimate for (4.3) follows by an adaptation of the spectral-calculus technique 
from [3, 4]. With this we prove the optimality of our approach to spectral asymptotic 
estimation. 

The following lemmata are the main technical results which are needed to estimate the 
quantities (4.3). 

Lemma 4.1. Take f eQZ "", then ^[H" 1 / - H^/] = (/, H" 1 /) - (/, H^/). 
Proof. The proof is a straight forward computation. Take / e " " , then 

h K [ulf} = (f,nlf) 

and we have 

h K [U-\f - Hlf] = (/, H" 1 /) - K{H-\f, Hlf) - h K (Hlf, H" 1 /) + (/, Hlf) 

= (/,h-V) - (h-vv.h^hL/) - (hV'h^/.h-vv) 
+ (/,hL/) 

= (/,H- 1 /)-(/,Ht /). 

□ 
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Lemma 4.2. Let f e Ti be given then set r/ :~ H fc / — Hj^H^/. //w assume 
||(Hy 2 H^ 1/2 )t|| < oo then q f = (H e /2 H- 1/2 )t r/ and 

t> b (KU,v) + (qf,n 1 e /2 v) = (f,v), v e Q. 
Furthermore, it holds that ||r/|| 2 = (/, H^ 1 /) - (^H^/). 

1 /2 

Prao/ It holds that r/ _L H 6 7 Qoo, which can be checked by a direct computation. The 
operator B := (H e ' H fe ) has the closed range so 

n - r(b*) e n(b) = r(b*) e h^q^. 

Therefore we have r/ € R(B*) and so we may write := BV/. A direct computation 
now shows that 

whu*) + fe.Hy 2 ,) = (hJ^h^/.h;/^) + (B*t r/I BH; /2 «) 

= (Hj/ 2 Ht / + r / ,Hj/ 2 t;) = (/,«). 

□ 

The main quantitative theorem about the asymptotic behavior of (4.1) follows now di- 
rectly. 

Theorem 4.3. Assume t := || (hV 2 Hj 1/,2 )t|| < oo then we have; for f <G Qoo ""; 

(4.4) (/.HrVM/.HL/) < (/)H _ 1/} _ (/;HL/) < ^((/.H^VH/.HL/)) 

and 

(4.5) -l^a.A^^tj^/cA 00 )^ ^7? 2 (0,A°°), i = l,...,m, 

where m is the multiplicity of the discrete eigenvalue A°° (not necessarily below the infi- 
mum of the essential spectrum ofTl^). 

Proof. For any / e Ti we have 

MHtj^ + te/iH^) = (/,«), neQ, 

h h (H- 1 /,«) + « 2 f) e (H- 1 /,t;) = (/,«), veQ. 

which implies 
and subsequently 

The right inequality in (4.4) follows from Lemma 4.2. To establish the left inequality 
of (4.4) we start from the identity [3, (22)]. We combine the integral representation for 
(/, H-\f) - (/, Hlf) from [3, pp. 41] and [3, (29)] to obtain 

(/,H- 1 /)-(/,Ht c /) = j o JT ^(dEn e (ml /2 ^ 1 f,K 1 e /2 Kb 1 f) 
= ((I + « 2 H e )- 1 H- 1 /,H- 1 /) 
> -^((I + He)- 1 ^- 1 /,^" 1 /) 

= -i((/,H r 1 /)-(/,HL/))- 
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The conclusion (4.5) for the approximation defects follows directly from the definition 
(3.16) and the observaton that 

holds. This completes the argument. □ 

Example 4.4. We will present this example as an abstract variation on (1.2). Let H be a 
positive definite operator, let P be a projection, R(P) C I?(H 1 / 2 ) and let := H^ 1 ^ 2 / — 
H^H^/. Consider 

t) K {u, v) = ((I + n 2 P)fi 1 / 2 u, a x ' 2 v) = t) b (u, v) + K 2 t) e (u, v), 

then 

ii(Hy 2 H- l/2 )t|| <i 

and (4.4) gives for / e N(PHe /2 ) 

r4fi , II^H 2 (f,H-if)-(f,H-V*P ± H-Wf) 1 

Ik/ll 2 (H-V2/,PH-V2/) - K 2- 

Here we have used || r f || 2 = (/, H," 1 /) - (/, H^/) and||r«|| 2 = (/, H" 1 /)-(/, H^/) 
to simplify the notation. On the other hand, we compute 

H- 1 =H- 1 / 2(P± + _L_ P)H -V2 

to establish 

(4.7) (/, H" 1 /) - (/, H-^P^H" 1 / 2 /) = ^—(H-V 2 /, ph- 1 / 2 /). 
Formulae (4.6) and (4.7) give 

1 _ Ik/ll 2 < J_ 



1 + k 2 \\r f \\ 

which is a very favorable estimate for k large. The lower estimate can be computed to be 



i lk^|| 2 



2 K 2 - ll^-ll 2 ' 

which is not as sharp as the upper estimate, but it is — newer the less — asymptotically 
optimal. 



5. A MODEL PROBLEM FROM ID THEORY OF ELASTICITY 

As an illustration of the applicability of Theorem 4.3, we consider the small frequency 
problem for the circular arch as described in [6, Chapter 8.8:3] and [28], cf. Figure 1. Let 
<t> : [0, /] — > R 2 be the middle curve of the arch. We take <j> to be the upper part of the 
circle with the radius R. The arch (the model problem we are considering) will be a thin 
homogeneous, elastic body of the constant cross-section A, whose area is A > 0. The arch 
will be clamped at one end and free at the other. The strain energy of the arch is given 5 by 



5 See also [34]. 
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FIGURE 1 . The Curved rod model 

the positive definite form 

' «<-) - e '[ (« + u i) (< (»; - i) (« - 1) *■ 

u,v £ Q{a) = {u£ H 1 [Q,l]x H 2 [0,l] : v(0) = 0, v' 2 (0) = 0}. 

Here u = (ui, v®) and v = (vi, V2) are the functions of the curvilinear abscissa s £ [0, 1], 
the constant E is the Young modulus of elasticity, the constant A is the area of the cross- 
section A and the constant / is the moment of inertia of the cross-section A . 

Let us assume we have the referent arch with the cross-section area A and the cross- 
section moment /. We consider the family of rods whose cross-section and the moment of 
inertia of the cross-section behave like 

A K = \a = e 2 A, I K = —1= s 4 I. 

We want to study the spectral properties of this family of arches as e — > 0. More general 
arch models can be treated by analogous procedures. This is a subject for future reports. 
For some given k > 0, k := £ _1 , we write 

- £ I'M + + l)' d » + ^ /(»'■ - 1) (« - 1) * 

and use A K to denote the operator which is defined by a K . Since A K has only the discrete 
spectrum we write A;(A K ), i £ N. After rescaling 

A,(A K ) = ~jK 

we see that Xf are the eigenvalues of the operator H K , which is defined by 
h K (u, v) = h 6 (u, v) + K 2 h e (u, v) 



RELATIVE CONVERGENCE ESTIMATES FOR THE LARGE COUPLING LIMIT 



21 



for u, v e Q(a K ) = <2(f) K ). Since Af enable us to describe only the eigenvalues of A K for 
which 

lim — A;(A K ) < oo. 

here we see where the name "low frequency problem", for the eigenvalue problem for H K , 
comes from. The low frequency problem satisfies the conditions of Theorem 2.1, so we 
conclude that the limiting form is 
(5.2) 

Wu,v) = EI J (u' 2 + ^-)' (v' 2 + ds, u, v e {f e Q(a K ),f[ - ^ = 0}. 

In [34] it has shown that (5.2) is the strain energy of the Curved Rod Model and that h K , 
k > are positive definite with 

Q(h K ) = {ueH 1 [0,l}xH 2 [0,l] : v(0) = 0, v' 2 (0) = 0}. 

Remark 5.1. From (5.2) we can see the significance of the condition 

(5.3) /( - ^ = 0. 

Assume the rod is locally straight. That is to say, assume R — > oo, then (5.3) turns into 

fx = 0, 

a condition of the inextensibility of the middle curve of the straight rod. The fact that 
f[ — ^ — is an inextensibility condition for the middle curve of the curved rod can be es- 
tablished by a rigorous differential geometric argument, see [34]. Continuing this heuristic 
reasoning, we conclude that Curved Rod model describes the transversal vibrations (per- 
pendicular to the middle curve) of the curved rod. Arch Model "couples" the longitudinal 
vibrations of the rod with the transversal vibrations. The study of finer properties of longi- 
tudinal vibrations requires the analysis of the so called "middle frequency problem", which 
will not be further considered here. However, since the "middle frequency problem" also 
falls under the scope of Theorem 2.1 this theory could also be applied in that case, too. 

5.1. Computational details. Based on (5.1) and (5.2) one concludes that the sequence 
h K satisfies the assumptions of Theorem 4.3. Here is a word of additional explanation 
in order. We have formulated all of our results about the forms t) b and h e based on the 
representations 

t> b (u,v) = (Hl /2 u,nl /2 v), 

However, we can represent (see (5.2)) the forms hf, and h e with the help of the operators 
R& : Q(f)b) — > "rib and R e : <2(h e ) — > TL e . The only assumptions on the operators 
(and R e ) is that they have a closed range in the auxiliary Hilbert spaces Hb (and H e ), cf. 
[16]. The representation theorem for the nonnegative definite forms implies 

(5.4) l) b (u,v) = (Hl ,2 u,iil ,2 v) = (R b u,R b v) nb , 

(5.5) *)e(u,v) = {H.l /2 u,U.l /2 v) = (R e ii,R e w)-H e , 

where (•, -)x generically denotes the scalar product in the Hilbert space X. The relations 
(5.4) and (5.5) imply that there exist isometric isomorphisms Q b : H b — > H and Q e : 



22 LUKA GRUBISIC 

1/2 1/2 

H e — > W such that H fe = QbRh, H e = Q e R e , and in particular 

(H^ /2 ?i, Hj /2 w) = (QfcR 6 u, QbRfeu) = (R&u, ~Rbv)n b , 
(Hy 2 u,Hy 2 w) = (Q e R e u,Q e R e w) = (R e u,R e v) We . 
We also have for u e <2(f)h) 

V2 _„ fWj(..> ^"iV 



Q^Hi^u = R fc u = VEl{u' 2 + 4 



Q^H^u = R e u = VEA (u\ ^) 

andR b : Q(h 6 ) -» W fc = L 2 [0, Z] and R e : Q(h e ) -» W e = L 2 [0, Z]. 

Note that H 6 is not positive definite but Hi, which is defined by the form hi = h h + h e , 
is. For the details see [19, 32]. If we were to change the notation we would have to set 
hfc := hi. Since this would unnecessarily complicate the exposition we opt not to do so. 

We show that 



I + AR 2 



(5-6) ||(Hy 2 H- 1/2 )t|| < y ( /?2 

for our model problem. Set t := || (He /2 H7 1/2 )t || then 

\mi /2 Ki 1/2 r«f\\= sup l( ^ 7 f : }l > \ \\p*~v\ 

vefl(W ||H/ v|| t 

since 



N((Hy 2 H7 1/2 )*) = N( H- 1/2 Hy 2 ) = N(H* /2 ) = Q x 

For Q- X q s e L 2 [0, Z] we define v = (/ ( ' ) (<5j 1 g/)(s)ds, 0) (an element of Q(h K )). For 
general v we have 



1/2 



Now, set v = v and compute 



1/2 VEI + EAR\ l 
" 1 V °" = 

This establishes 

|( g/ ,H* /2 v)| |( g/ ,H* /2 v )| KQJ^R^VqM _ / Ai? 2 
v5S>) ||H;/ 2 v|| " ||Hi/ 2 v || - v^CT ||g/|| Vi + Ai? 2 ^ 11 ' 
which completes the proof of (5.6). 

5.2. Quantitative (and qualitative) conclusions. The fact (5.6) allows us to apply The- 
orem 4.3 to obtain precise estimates for the behavior 6 of rji(e, A£°). Since Hi and not 
H;, is the positive definite operator, we will use the rod with the diameter e as a referent 
configuration. We chose 



( j ° _ 6 V A& A-+A-' 

where A^ c and Aj^ n are the two lowermost eigenvalues of the Curved Rod model and 
A^ in denotes the lowermost eigenvalue of the Arch Rod Model of the rod with diameter 



6 We have tacitly dropped the exponent from r\i (e 2 , A2°) in order to simplify the notation. 
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e. Theorems 3.11, 4.3, Remark 3.6 and Corollary 3.10 — together with the observation that 
r]i(e, A°°) < 1 for any e > — directly imply that 

~ ^min ^ _2 + A R 2 ) Agg c + A^j n 

Too ^ £ Tm xoo _ X oo > 0<e<e O - 

A min A 11 A sec A m in 



Furthermore, if we chose e\ = if \J 'Wi - Ag| + a~ » then we obtain 

2 2(J + ^ i? 2 )r/,( £l , AgjJ A- - A £ min 4(1 + A R 2 ) A s °g c + A~ 
( j ^ 2 " A- " -4 R? \? ec -\% n > 

< £ < E\. If we are only interested in the upper estimate and we assume that there is 
m e N such that A^ < A~ +1 , then we have 

(59 ) V°- X ^ ? MI + AR 2 ) % _ 

{ ^ } Af S . |A^-A-| AR 2 £ ' > • • • > 

1 max mm — — — 

i=l,...,m kjti A^ + A°° 

This estimate holds for all e < £2, where 62 is defined as the first e for which the righthand 
side of (5.9) is smaller then 1. Estimate (5.9) can naturally be refined with the use of other 
Vi( £ 2, A°°) as is given by the framework of Theorem 3.3. The optimality of the estimate is 
meant in the sense of Theorem 3.9. 



6. Conclusion 

We have presented a constructive approach to spectral asymptotic estimates in the large 
coupling limit. Although we have concentrated on a use of this results for theoretical 
considerations from [5, 7, 9, 28], they are expected to be particularly useful in a design of 
computational procedures for various singularly perturbed spectral problems. This can be 
illustrated when comparing the numerical procedures for the Arch Model and the Curved 
Rod Model. It has been shown that the Curved Rod Model is is better behaved, with 
respect to the finite element approximations than the Arch Model, see [33]. Furthermore, 
a qualitative conclusion of Section 5.2 is that when interested in the transversal vibrations 
only, Arch Model can be ignored (up to the corrections of order e 2 ). For more on the lower 
dimensional approximations in the theory of elasticity see [6, 19, 28, 32, 34]. 

In a practical computational setting it is not reasonable to assume that the spectral prob- 
lem for Hqo will be exactly solvable. We would like to emphasize that in the design of this 
theory we have not built the requirement of the explicit solvability of Hoc into our results. 
To be more precise, nowhere in the proofs of Theorems 3.5 and 3.11 or Corollaries 3.7 and 
3.8 is it necessary to have R(P) = £oo- The only place where this assumption was neces- 
sary was to establish that (3.6) and (3.16) define the same approximation defects. Theorem 
3.3 and similar results from [13, 17] — which are the workhorses of this theory — do not 
need this assumptions. Subsequently, the only limiting factor is the computability of rji(P) 
and the availability of information on the distance of spec(S K ) — from Theorems 3.2 and 
3.3 — to the unwanted component of the spectrum. With this we hope to have illustrated 
the advantages and limitations of our theory 
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